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Supersaturated design (SSD) has received much recent interest 
because of its potential in factor screening experiments. In this pa- 
per, we provide equivalent conditions for two columns to be fully 
aliased and consequently propose methods for constructing -B(/nod)- 
and x^-optimal mixed-level SSDs without fully aliased columns, via 
equidistant designs and difference matrices. The methods can be eas- 
ily performed and many new optimal mixed-level SSDs have been ob- 
tained. Furthermore, it is proved that the nonorthogonality between 
columns of the resulting design is well controlled by the source de- 
signs. A rather complete list of newly generated optimal mixed-level 
SSDs are tabulated for practical use. 

1. Introduction. The supersaturated design (SSD) is a factorial design 
in which the number of runs is not sufficient to estimate all the main effects. 
Such designs are useful when the experiment is expensive, the number of 
factors is large, and only a few significant factors need to be identified in 
a relatively small number of experimental runs. Booth and Cox (1962) first 
examined these designs systematically and proposed the E{s'^) criterion. 
However, such designs were not further studied until the appearance of the 
work by Lin (1993, 1995), Wu (1993), Tang and Wu (1997) and Cheng 
and Tang (2001). Research on mixed- level SSDs includes the early work by 
Fang, Lin and Liu (2000, 2003) who proposed the -E'(/nod) criterion and the 
FSOA method for constructing mixed-level SSDs, and work by Yamada and 
Matsui (2002) and Yamada and Lin (2002) who used to evaluate mixed- 
level SSDs. Recent work on mixed-level SSDs includes Xu (2003), Fang et al. 
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(2004a), Li, Liu and Zhang (2004), Xu and Wu (2005), Koukouvinos and 
Mantas (2005), Liu, Fang and Hickernell (2006), Yamada et al. (2006), Ai, 
Fang and He (2007), Tang et al. (2007), Chen and Liu (2008a, 2008b), Liu 
and Lin (2009), Liu and Cai (2009) and Liu and Zhang (2009). 

This paper proposes some methods for constructing i?(/NOD)- and 
optimal mixed-level SSDs without fully aliased columns, and with a control 
on the nonorthogonality. A large number of optimal designs is obtained. The 
remainder of this paper is organized as follows. Section 2 provides relevant 
notation and definitions. In Section 3, we propose the general construction 
methods for mixed-level SSDs along with illustrative examples. Discussions 
on the nonorthogonality of the resulting designs are given in Section 4. In 
Section 5, a review of the existing methods for mixed- level SSDs and compar- 
isons with the current methods are made, and some concluding remarks are 
provided. For coherence of presentation, all proofs are placed in Appendix A 
and newly constructed designs are tabulated in Appendix B. 

2. Preliminaries. A mixed-level design that has n runs and m factors 
with qi,...,qm levels, respectively, is denoted by F{n,qi ■ ■ ■ Qm)- When 
Yl^=iilj ~ 1) = ~ li the design is called a saturated design, and when 
Z]j=i(9j — 1) > "ra — 1, the design is called a supersaturated design (SSD). 
An F(n,qi ■ ■ - qm) can be expressed as an n x m matrix F = (fij). When 
some qj^s are equal, we use the notation F{n,q[^ ■ ■ ■ q^') indicating rj fac- 
tors having qi levels, i = 1, . . . ,1. If all the g^-'s are equal, the design is said 
to be symmetrical and denoted by F{n,q"^). Let /j be the ith row of an 
F{n, qi - ■ ■ qm) and be the jth column which takes values from a set of qj 
symbols {0, . . . ,qj — 1}. If each column is balanced, that is, it contains the 
qj symbols equally often, then we say F is a balanced design. Throughout 
this paper, we only consider balanced designs. Two columns are called fully 
aliased if one column can be obtained from the other by permuting levels; 
and called orthogonal if all possible level-combinations for these two columns 
appear equal number of times. An F(n, qi - ■ ■ qm) is called an orthogonal ar- 
ray of strength two, denoted by Ln{q^) for the symmetrical case, if all pairs 
of columns of this design are orthogonal. 

The set of residues modulo a prime number p, {0, 1, . . . ,p — 1}, forms a field 
of p elements under addition and multiplication modulo p, which is called 
a Galois field and denoted by GF(p). Note that the order of a Galois field 
must be a prime power. A Galois field of order g = for any prime p and any 
positive integer u can be obtained as follows. Let g{x) = 6o + + • ■ • + h^x"^ 
be an irreducible polynomial of degree u, where bj € GF(p) and = 1. 
Then the set of all polynomials of degree n — 1 or lower, {gq + aix + • • • + 
au-ix'^~^\aj S GF(p)}, is a Galois field GF(g) of order q = p'^ under addition 
and multiplication of polynomials modulo g{x). For any polynomial f{x) 
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with coefficients from GF(p), there exist unique polynomials q{x) and r{x) 
such that f{x) = q{x)g{x) + r{x), where the degree of r{x) is lower than u. 
This r(x) is the residue of /(x) modulo g{x), which is usually written as 
f{x) = r{x){modg{x)). 

A difference matrix, denoted by D{rq, c, q), is an rqx c array with entries 
from a finite Abelian group {A, +) with q elements such that each element of 
A appears equally often in the vector of difference between any two columns 
of the array [Bose and Bush (1952)]. Note that if A is an Lrq{q'^), then 
it is also a difference matrix. A difference matrix D{rq,c,q) with c > 1 is 
said to be normalized, denoted by ND(rg,c, g), if its first column consists 
of all zeros. In fact, for any difference matrix D, if we subtract the first 
column from any column, then we can obtain a normalized difference ma- 
trix. 

For a scalar a and a matrix A, let a + A denote the element- wise sum of 
a and A. For any two matrices A = (aij) of order r x s and B of order ux v, 
their Kronecker sum and Kronecker product are defined to be 

(aii+B ■■■ ais + B\ / auB ■■■ auB' 
and A^B 
Clfi -\- B ■ ■ ■ drs ~l~ / \ CLrlB ■ ■ ■ d^sB 

respectively. Here, we use "+^" and "©^" to denote the sum and Kronecker 
sum defined on A, respectively. 
For a design F = {fij)nxm, let 

m 

A,,(F) = j;5|f and u^,,{F) = J2QkS\-\ 

k=l k=l 



where 6^j^ = 1 if /j^ = /j^, and otherwise. Then Xij{F) and Uij{F) are 
called the coincidence number and natural weighted coincidence number be- 
tween rows fi and fj , respectively. A design with equal coincidence numbers 
between different rows is called an equidistant design. From Mukerjee and 
Wu (1995), a saturated Ln{q^) is an equidistant design with 

TTl — 1 

(1) Xij[F) = and uJij{F) =m—l for i ^ j. 

The -E'(/nod) criterion proposed by Fang, Lin and Liu (2000, 2003) is 
defined to minimize 

2 ^ 

-E(/nOd) = 7 7T 2^ fNOD{f\f), 

mlm — 1 ^-^ 

^ ' l<i<j<m 



where 



n 



/NOD(r,/^) = ^5;(n„fe(r,/ 

a=0 6=0 ^ 
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nab{f '',f'^) is the number of (a, 6)-pairs in (/*,/-'), and n/{qiqj) stands for 
the average frequency of level-combinations in (/*,/-'). Here, the subscript 
"NOD" stands for nonorthogonality of the design. The /nodC/*,/-') value 
gives a nonorthogonality measure for (/*,/-'), and columns /* and are 
orthogonal if and only if /nod(/*)/'') = 0. It is obvious that F is an or- 
thogonal array if and only if -E'(/nod) = 0, that is, /nodI/*,/-') = for all 
i,j = 1, . . . ,m,i ^ j . Thus -^(/nod) measures the average nonorthogonality 
among the columns of F. 

Another criterion that is to be minimized was defined by Yamada and Lin 
(1999) and Yamada and Matsui (2002) as x^{F) = Y.i<i<j<m1i1jhoD{f, 
f^)/n. Obviously, ^'(/nod) and x^i^) are equivalent in the symmetrical 
case. Here, we adopt both -^'(/nod) and x^{P) to evaluate the newly con- 
structed SSDs. There are also some other criteria for assessing mixed-level 
SSDs [see, e.g., Liu and Lin (2009) for a general review]. 

The following results, regarding the -^(/nod) and x^{F) optimality cri- 
teria of a design, will be needed for our construction methods. 

Lemma 1. (a) [Fang et al. (2004a)]. If the difference among all coinci- 
dence numbers between different rows of design F does not exceed one, then 
F is E{fjsioi)) -optimal. 

(b) [Li, Liu and Zhang (2004); Liu, Fang and Hickernell (2006)]. // the 
natural weighted coincidence numbers between different rows of design F 
take at most two nearest values, then F is -optimal. 

3. Proposed construction methods. In this section, we first provide some 
equivalent conditions for two columns to be fully aliased, then propose meth- 
ods for constructing i?(/NOD)- and x^-optimal SSDs, and finally study the 
properties of the resulting designs. 

3.1. Equivalent conditions for two columns to be fully aliased. An 
£J(/nod)- or x^-optimal SSD may contain fully aliased columns, which is 
undesirable. Let matrix Xj = (x^J of order n x qj be the induced ma- 
trix [Fang et al. (2004a)] of the jth column of an F{n,qi ■ ■ ■ qm), that is, 
= 1 if the sth element in the jth column is t — 1, otherwise 0, for 
s = 1, . . . ,n,t = 1, . . . ,qj and j = 1, . . . ,m. The following theorem presents 
theoretical results concerning the column aliasing that will be used in the 
construction methods. 

Theorem 1. Suppose Xj = (x^^) is the induced matrix of a balanced 
column f^ = {fij, . . ■ , fnjj)' with qj levels, j = 1,...,4, and rii = n^, 

712 = "-4- 
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(a) For qi=q2=qs = qi = q and A = {0, . . . ,q - I}: 

(i) and f ^ are fully aliased if and only if XiX[ = X^X!^; 

(ii) the induced matrix of is . . . , (X2P/„^J']' = 

{Xi (g) X2)P, where P = {Pq, . . . ,Pg_iy and Pi is a permutation matrix 
defined by 

i+A{0,...,q-l) = {0,...,q- l)Pl, i = 0, ...,(?- 1; 

(iii) iff^^Af"^ o.nd f^ (B^f^ are fully aliased, then f^ is fully aliased 
with f^ and is fully aliased with 

(b) (i) The induced matrix of the qiq2-level column q2{f^ — ^^Y^)(B if — 

(ii) columns q2{f^ - ^) © (/^ - 32^) + and q^f^ - © 
_ 94-1 ^ _|_ 9394-1 ^j.g /ti/Zy aliased if and only if f^ is fully aliased 

with f^ and is fully aliased with 

(iii) for q3=q, = q, <?2(/^ -'-^)®{f^-'-^) + and f^ ®A f 
are not fully aliased in any case. 

3.2. Construction of optimal symmetrical SSDs. We next present the me- 
thods for constructing -^(/nod)- and x^-optimal SSDs without fuhy aliased 
columns. 

Theorem 2. Let D be an ND(rg,c, (7) defined on an Abelian group A = 
{0, . . . ,q — 1} without identical rows, F be an F(n, q^) without fully aliased 
columns and with constant coincidence numbers, say X, between its different 
rows, then: 

(a) F(BaD' is an F{cn,q^^'"^) with two different values of coincidence 
numbers, mr and Xrq; 

(b) F (Ba D' has no fully aliased columns. 

From Lemma 1, if \mr — Xrq\ < 1, then F (Ba D' is both -E(/nod)- and 
X^-optimal. The following corollary can be directly obtained from Lemma 1, 
Theorem 2, and equation (1). 

Corollary 1. Let F be a saturated Ln{q^) and D be an ND(g,c, g) 
without identical rows. Then F (Ba D' is an F{cn, q^'^) without fully aliased 
columns and with two different values of coincidence numbers, m and m — 1, 
and thus is both -^(/nod)- and -optimal. 
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Table 1 
F and D in Example 1 
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From Hedayat, Slone and Stufken (1999), there exist an Ln{q"^) with n = 
and m = (n — 1) /(g — 1) and an ND(g, g, q) without identical rows for any 
prime power thus optimal F{cq^ ^q^'^^^^~^'^/'^'^~^^) designs with coincidence 
numbers (g* — — 1) — 1 or (g* — l)/{q — 1) can be constructed from 

Corollary 1, where c is a positive integer and c< q. 

Example 1. Let F be an ^9(3"^) and D be an ND(3,2,3) (cf. Table 1), 
then F(BaD' is an F(18,3^^) with coincidence numbers 4 and 3 as listed in 
Table 2, where A = GF(3). This new design is an -E'(/nod)- and x^-optimal 
SSD without fully aliased columns. 
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The _F(18, S"'^'^) constructed in Example 1 
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Table 3 
Fi , F2 and D in Example 2 
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3.3. Construction of optimal SSDs with two different level sizes. Based 
on Lemma 1 and Theorem 2, the following theorem can be obtained. 

Theorem 3. Let Fi he an F{ni,q^^) with constant coincidence numbers 
Xi, and no full aliased columns, i = 1,2. Let D be an ND(rgi, n2, defined 
on Abelian group Ai = {0, — 1} without identical rows. Then F = 

(-^1 ffi^i L)',Oni ® F2) is an F{nin2,qi"^^'^^q^^) without full aliased columns. 
Furthermore: 

(a) if |(A2 + rmi) — (m2 + Airgi)| < 1, then F is E{f^oii)-optimal; 

(b) «/(?2A2 + qirmi = (?2"i2 + Xirq\, then F is -optimal. 

Next, let us consider two illustrative examples for Theorem 3. 

Example 2. Let Fi be an Li{2^), F2 be the £'(/NOD)-optimal F(6,3^) 
obtained by Fang, Ge and Liu (2004) and D be an ND(8,6, 2) without 
identical rows obtained from an Ls{2'^) based on A = GF(2). Then Ai = A2 = 
l,qi = 2,g2 = 3,r = 4, mi = 3 and 7712 = 5 which satisfy the condition that 
A2 + rmi = 7712 + Air(7i = 13, thus (Fi ©_4 D', O4 © F2) is an F(/NOD)-optimal 
F(24, 2^^3^) with constant coincidence numbers 13. The source designs and 
resulting design are listed in Tables 3 and 4, respectively. 

Example 3. Let Fi be an ^4(2^), F2 be the F(6,3^°) obtained by 
Georgiou and Koukouvinos (2006) and D be an ND(24, 6, 2) without iden- 
tical rows obtained from an L24(2^^) based on ^ = GF(2). Then Ai = 1, 
A2 = 2,gi = 2,52 = 3,r = 12,7771 = 3 and m2 = 10 which satisfy the con- 
dition that q2X2 + (71 77771 = '?2"72 + ^i^qf = 78, thus {Fi (Ba -D',04 © F2) 
is a x^(-C))-opti™cil ^(24,2^^3^^*) with constant natural weighted coinci- 
dence numbers 78. Exact details are omitted here but available upon re- 
quest. 
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Table 4 

The _F(24,2^''3^) constructed m Example 2 
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3.4. Construction of optimal SSDs with three different level sizes. The 
next lemma is useful in the upcoming proposed construction method. 

Lemma 2. Let F = {-2^, -2^, . . . , 2^, 2^} = {0, . . . , g - 1} - 2^ 

and Vi = {i - + V,i = 0, . . . ,p - 1, then ViDVj = ^ for i ^ j and 

Ur=o^^ = {-^'-^--->^>^} = {0,---,P9-l}-2^, where'^ 
is an empty set. 

From this lemma, we can obtain the following theorem in a straightfor- 
ward manner. 

Theorem 4. Let Fi be an F{ni,q^'^) with constant coincidence numbers 
Ai, i = 1, 2, then 92(^1 - ^)® {F2 - + %s an F{nin2, {qmT''''^) 

with three different values of coincidence numbers Xim2, X2'itii and A1A2. 
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This theorem, along with Lemma 1 and Theorem 2, leads to the following 
theorem, which provides another construction method of -E(/nod)- and 
optimal SSDs. 

Theorem 5. Suppose Fi is an F{ni,q^^) with constant coincidence 
numbers Aj and no fully aliased columns, i = 1, . . . ,4, is an ND(r3g3, n2, qs) 
defined on Abelian group ^3 = {0, . . . , — 1} without identical rows, D4 is 
an ND(r4g4, ni, 94) defined on = {0, . . . ,q4 — 1} without identical rows, 
and they satisfy (i) ui = n^,n2 = n^; (ii) the first rows of and F4 consist 
of all zeros; (iii) there are no fully aliased columns between F3 and D'^ or 
between F4 and D'^ . Then 



is an F{nin2,{qiq2)"^^"^^Q3^^^^'^'^Q4^^^'^'^*) without fully aliased columns and: 

(a) if the difference among three values X2mi + r3m3 + A4r4(74,AiA2 + 
'''3^3 + ?'4?Ti4 and \im2 + ^sr^qs + r4m4 does not exceed one, then F is 
E{f ^od) -optimal; 

(b) if qiq2hmi + 93?'3"^3 + A4r4g| = g'iq'2AiA2 + q3r3m3 + g'4r4m4 = 
qiq2^i'm2 + \3rzq^ + qir^m^^, then F is -optimal. 

The following two examples serve as illustrations of the construction 
method in Theorem 5. 

Example 4. Let Fi and F3 be two L4{2^ys; F2 be the ^(6,2^°) ob- 
tained by Liu and Zhang (2000); F4 be the F(6,3^) obtained by Fang, 
Ge and Liu (2004); D3 be an ND(12,6,2) without identical rows obtained 
from an Li2i2^^); D4 be an ND(12,4,3) without identical rows; ^3 = GF(2) 
and A4 = GF(3). Suppose the first rows of F3 and F4 consist of all ze- 
ros. Then based on Theorem 5, Ai = A3 = A4 = 1, A2 = 4, mi = = 3, m2 = 
10, 777,4 = 5, g'l = 92 = 93 = 2, 94 = 3, r3 = 6, r4 = 4 and A2mi -|- r3m3 -|- A4r4g4 = 
AiA2 + r3m3-|-r4m4 = Aim2 + A3r3g3-|-r4m4 = 42. Thus, from (2), we obtain 
an £'(/NOD)-optimal F(24, 4^'^2^^3^'^) with constant coincidence numbers 42 
and no fully aliased columns. 

Example 5. Let Fi and F3 be two L4(2^)'s; both F2 and F4 be the 
F(6, 3^) obtained by Fang, Ge and Liu (2004); D3 be an ND(24, 6, 2) without 
identical rows obtained from an L24(2^^) based on A3 = GF(2) and D4 be an 
ND(6,4,3) without identical rows based on A4 = GF(3). Suppose the first 



(2) 




^^3 0^3 D'3,D'^ ®A,Fa 
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rows of F3 and i<4 consist of all zeros. Then Ai = A2 = A3 = A4 = 1, = 93 = 
2, q2 = qA = 3, = 12, = 2, mi = = 3, m2 = = 5, which satisfy the 
condition that qiq2\2'mi + qsrsm^ + X^r^ql = qiq2XiM + q'ir^m^ + qAr^mi = 
(7i(72Aim2 + A3r3^| + q4r4m4 = 108. Thus, the design constructed through (2) 
is a x^-optimal F(24, 6^^2^^3^'^) with constant natural weighted coincidence 
numbers 108 and no fully aliased columns. 

4. Nonorthogonality of the resulting designs. In the previous section, 
construction methods for -E(/nod)- as well as x^-optimal SSDs without fully 
aliased columns are provided. Full aliasing can be viewed as the extreme case 
of nonorthogonality. In this section, we will investigate nonorthogonality, 
measured by /nod, of the resulting designs, and show how it is controlled 
by the source designs. 

Theorem 6. Suppose /* = {fu, ■ ■ ■ , fmi)' is a qi-level balanced column 
with induced matrix Xi, Ai = {Q, . . . ,qi — 1}, i = 1, . . . , 4, rii = 723, 712 = 714. 
Let h, = q2if' -^)(B (f -^) + and /12 = q^if - (/^ - 

2^) + 2^. Then: 

Ja) hoDihl,h2) = fNOD{f\f)fNOD{f,f) + Jj/nOd(/\/') + 

^/nod(/^,/^); 

(b) if qi = q2, 93 = 94, then 

/nOd(/' f,f ®As h < qiq3fNOD{f\f)fNOD{f,f) 

+ mm{nlhoB{f\f),nlhoB{f, /')}, 

where the equality holds if and only if f^ is orthogonal to f^ or is orthog- 
onal to f^; 

(c) ifqi = q2, then 

/nod(/' ®A, f\h2) < 9i/nod(/\/')/nod(/',/') 

+ min{n2/94/NOD(/\ /^), ni/93/NOD(/^, /^)}, 

where the equality holds if and only if f^ is orthogonal to f ^ or f^ is orthog- 
onal to 

Theorem 6 shows that the nonorthogonality measured by /nod of the 
resulting designs is well controlled by the source designs. If the source designs 
have small values of /nod, then the resulting design will also have small 
values of /nod- In particular, we have the following. 

Corollary 2. Suppose /'* = {fu, . . . , fn^i)' is a qi-level balanced column 
with induced matrix Xi, Ai = {Q, . . . , qi — 1} , i = 1, . . . ,4, ni = 723,712 = 714. 
Then: 
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(a) if is orthogonal to and is orthogonal to , then 92 (/"^ — 
e (/2 - 22^) + IS orthogonal to q^if - 22^) (/^ - 3^) + 

9394-1 . 
2 ' 

(b) ifqi = q2, Q3 = Qi o-nd is orthogonal to or is orthogonal to f , 
then (Bai orthogonal to ©^3 

(c) if qi = q2 and f^ is orthogonal to f^ or f^ is orthogonal to f^, then 
f^ ®Ai f^ is orthogonal to q^if^ - (/^ - + 2^^. 

This corollary indicates that the orthogonality between columns of the 
source design is maintained in the generated designs. 

5. Discussion and concluding remarks. In this paper, we have presented 
some construction methods for i?(/NOD)- and x^-optimal SSDs. A review of 
the existing methods for mixed-level SSDs and comparisons with the current 
methods are summarized below. 

(a) Yamada and Matsui (2002) and Yamada and Lin (2002) proposed two 
methods for constructing mixed- level SSDs consisting of only two- and three- 
level columns through computer searches. However, their resulting designs 
have no theoretical support and typically are unable to achieve the lower 
bound of x^-value. 

(b) Fang, Lin and Liu (2000, 2003) proposed an FSOA method for con- 
structing £^(/NOD)-optimal mixed- level SSDs from saturated orthogonal ar- 
rays. Li, Liu and Zhang (2004) and Ai, Fang and He (2007) extended the 
FSOA method to construct x^-optimal SSDs. Koukouvinos and Mantas 
(2005) constructed some £'(/NOD)-optimal mixed-level SSDs by juxtapos- 
ing either a saturated two-level orthogonal array and an £^(/NOD)-optimal 
mixed-level SSD, or two E(/NOD)-optimal SSDs. Fang et al. (2004a) and 
Tang et al. (2007) presented some methods for constructing -B(/nod)- and 
X^-optimal mixed-level SSDs, respectively, from given combinatorial designs. 
There are many constraints on the parameters of saturated orthogonal ar- 
rays and combinatorial designs and the construction of most combinatorial 
designs are unresolved. Thus, the optimal SSDs obtained by their methods 
are rather limited. 

(c) Yamada et al. (2006) presented a method for constructing mixed- 
level SSDs by juxtaposing two SSDs, each of which is generated by the 
operation "0" of an initial matrix and a generating matrix. It can be seen 
that their operation "©" is in fact equivalent to the "0^^" in this paper with 
A= {0,. . . ,q — l}, and they only provided the theoretical justification of the 
X^-optimality for the SSD with n = 6. Recently, Liu and Lin (2009) proposed 
a method to construct x^-optimal mixed-level SSDs from smaller multi- level 
SSDs and transposed orthogonal arrays based on Kronecker sums. It can be 
easily confirmed that the result of Liu and Lin (2009) is merely a special case 
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of our Theorem 3, by taking Fi as Lg-^ = (0, . . . , gi — 1)' and D as L^gi (^i^)- 
Thus, ah their designs can be constructed by our Theorem 3. 

(d) Using fc-cychc generators, Chen and Liu (2008a) and Liu and Zhang 
(2009) constructed some -E'(/nod)- and x^-optimal mixed-level SSDs, re- 
spectively. The /c-cyclic generators were obtained via computer searches, 
when the values of k, the run size and/or the level sizes become larger, the 
computer searches tend to be ineffective and impractical. 

(e) Recently, Liu and Cai (2009) proposed a new construction method, 
called the substitution method, for E'(/NOD)-optimal SSDs. It can be seen 
that all the i?(/NOD)-optimal SSDs tabulated in our Tables 6 and 8 are 
different from those tabulated in their Appendices. 

Note that the newly proposed methods use small equidistant designs and 
difference matrices to generate large designs. Many difference matrices can 
be found in Hedayat, Slone and Stufken (1999), Wu and Hamada (2000) and 
from the site http:/ /support. sas.com/techsup/technote/ts723. html maintai- 
ned by Dr. W. F. Kuhfeld of SAS. Equidistant designs can be found in 
Ngugen (1996), Tang and Wu (1997), Liu and Zhang (2000), Lu et al. 
(2002), Fang, Lin and Liu (2003), Fang, Ge and Liu (2002a, 2002b, 2004), 
Lu, Hu and Zheng (2003), Fang et al. (2003, 2004a, 2004b), Aggarwal and 
Gupta (2004), Eskridge et al. (2004), Georgiou and Koukouvinos (2006), 
Georgiou, Koukouvinos and Manias (2006), Chen and Liu (2008a), Liu and 
Cai (2009) and others. Difference matrices can also be obtained from or- 
thogonal arrays or by taking the Kronecker sums of difference matrices, 
and equidistant designs also include saturated orthogonal arrays of strength 
two. 

The appealing feature of our methods is that they can be easily applied 
and the resulting designs are £'(/nod)- and/or x^-optimal SSDs without 
fully aliased columns. In particular, the nonorthogonality between columns 
of the resulting designs is well-controlled by the source designs, that is, if 
the source designs have little nonorthogonality, the generated design will also 
have little nonorthogonality. From these proposed methods, many optimal 
SSDs can be constructed in addition to those tabulated in the Appendix. 

In regard to the statistical data analysis for mixed-level SSDs, it should 
be noted that analyzing the data collected by such SSDs is a very important 
but complicated task which has attracted much recent attention. See, for ex- 
ample, Zhang, Zhang and Liu (2007), Phoa, Pan and Xu (2009) and Li, Zhao 
and Zhang (2010). When there are many more factors than the number of 
permitted runs due to expense (e.g., money or time), the nonorthogonality 
among factors may be very severe and may prevent the few active factors 
to be identified correctly by any existing method. Therefore, the data anal- 
ysis for SSDs in general remains an important and challenging topic for 
further research. Some recent study on the analysis of "high-dimension and 



CONSTRUCTION OF MIXED-LEVEL SUPERSATURATED DESIGNS 



13 



low-sample size" in genetic studies (e.g., studying 6,000 genes with only 37 
observations) may be relevant. 



APPENDIX A: PROOFS 

Proof of Theorem 1. (a)(i) If and f '^ are fully aliased, that is, 
can be obtained by permutating the levels of then there must exist 

a permutation matrix Q of order q that satisfies Xi = X^Q, thus XiX[ = 

XsQQ'X'^ = X^X'r^. 

On the other hand, let Vi and V3 be the vector spaces spanned by the 

columns of Xi and X3, respectively. If XiX[ = X^X^, then 14 = V3, and for 

any column of Xi , we have 



x° = kixl + 



where x\ is the ith column of X3 , i 



l,...,q. 



Since any two columns in an induced matrix share no element 1 at any 
position and each column has ni/q ones and ni — ni/q zeros, there must 
exist only one ki 7^ 0, that is, x^ is identical to a column of X3. Then there 
exists a permutation matrix Q of order q satisfying Xi = X3Q, thus can 
be obtained by permutating the levels of f^, that is, and are fully 
aliased. 

(h) Note that the induced matrix of f is )',..., (X2P/„^ J']' 



andP/,,=E?= 



^xl^Pt-i, then 



X9P 



hi 



:diag{X2,...,X2} 



^X2Pf„^l. 



n-i 



t=l 



{In, ^ X2)iXi Ig 



{Xi 



i-1 



t=l 
P = 



X2)P. 



where /„ is the identity matrix of order n. 

(hi) The induced matrices of and ©^ are {Xi 0X2) P and 

(X3 (g) X4)P, respectively. If ©^ and ©^ are fully aliased, then 
there exists a permutation matrix Q of order q such that 

{Xi ® X2)P = (X3 A:4)Pg, that is, Y,^ltX2Pt-i = ^xltX4Pt-iQ 



t=i 



for s = 1, . . . ,ni. For any s and i, there is only one nonzero element of x*^ 
for t = 1, . . . ,g that equals 1. Thus, X2Pt^_i = X^Pt^^iQ, and therefore 
is fully ahased with Similarly since ©^ and ©^ are also fully 
aliased, it follows that is fully aliased with f ^. 
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(b)(i) It can be obtained easily from the definition of an induced matrix. 

(ii) Prom (a)(i), we only need to prove that columns q2{f^ - ^) (f - 
22_d) + 21^ and q4{f-^)(B{f^-^) + are fully ahased 
if and only if XiX[ = X^X'.^ and XsX^ = X^X'^. From (b)(i), the induced 
matrices of these two columns are Xi ® X2 and ® X4 , respectively, thus 
from (a)(i), they are fully aliased if and only if {Xi ®X2){Xi (8X2)' = {X^ (8) 
X4){X3 (g) Xi)', that is, XiX{ (g) X2X2 = X3X3 (g) which means that 
XiX[ = 0X3X3 and X2X2 = 1/0X4X4 for some a 7^ 0. Since the elements in 
XiX- are all ones and zeros for i = 1, . . . , 4, then a = 1, that is, XiX( = X3X3 
and XaX^ = X4X^. 

(iii) The induced matrices of columns q2{f^ — ^^^^) © (/^ — ^^^) + '^^'^2~^ 
and f ®A f are Xi X2 and [{XiPf^J , {XiPf^_^J]', respectively. If 
these two columns are fully aliased, then ^4^X2X2 = X4Pf.^Pj_^X'^, where 
Vij is the (i,j)th entry of XiX{,i,j = l,...,ni. Note that Vij can be zero, 
and hence ^4^X2X2 can be a zero matrix which contradicts the fact that 
X4^Pf..^P'jr,^X'^ cannot be a zero matrix in any case. □ 

Proof of Theorem 2. (a) Consider the ith and jth rows of F (Ba 
D', (/ii ®A di^y and ®a dj^)' , where i = {h - l)c + 12, j = (j'l - l)c + 
j2,h,ji = 1, • • • ) ^2,^2 = 1, • • • , c, and i 7^ j, fk and (i^ are the kth rows of 
F and D' , respectively. Then the coincidence number between (/j^ ©71^12)' 
and (/j^ (Ba '^j^)' equals the number of zeros in (/j^ — /jj (Ba {di^ — dj^). 

(i) Suppose ii = ji,«2 7^ j2i then /jj = /j^ and di^ ¥'dj2- From the defi- 
nition of difference matrix, each element in A occurs r times in di^ — dj^ . 
Therefore, (/jj — /jj) (Ba {di2 — dj2) = Om ®a {di2 — dj^) and there are mr 
zeros in Om ®a {di2 — dj2 ) , where 0^ denotes the m x 1 column vector with 
all elements zero, that is, the coincidence number between (/j^ ©^^^3)' and 
ifji ®Adj2y is mr. 

(ii) If ii ^ ji,i2 7^i2, similar to (i), it can also be easily seen that there 
are mr zeros in (/jj — /jj (djj — (ijj), that is, the coincidence number 
between (/jj ©^djj)' and {fj^ ©y^djj)' is mr. 

(iii) If ii / ji, Z2 = j2, that is, fi^ ^ fj^ , = dj^, then (/j^ - /j J ©^ (dij - 
•^j^) = ifii ~ fji) ®A Org, and there are Xrq zeros in (/j^ — /jj (Ba ^rg, that 
is, the coincidence number between (/j^ ©yidjj)' and [fj^ ©y^^j^)' '^''Q- 

(b) F (Ba D' can be obtained from D' (Ba F through row and column 
permutations. Thus, if D' (Ba P has no fully aliased columns, neither does 
F ®aD'. Let ©^ and ®a P be two different columns of D' ©^ F, 
where = {du, . . . , dd)' and /* = {fu, . . . , /„j)' for z = 1 and 2 are columns 
of D' and F, respectively. Since D is a normalized difference matrix, du = 
for i = 1 and 2. Let Xi and X2 be the induced matrices of and 
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respectively. Then the induced matrices of (Ba and (Ba P a^re 



Suppose 
(3) 



A and 



/ ^2 \ 



and 



V 



respectively. 



^A are fully aliased. Then from Theorem 1, 



Noting that Pd^^ = Pd^^ 
f = 1 in (3): 



^iPduP'dji^l - ^2PdaPdj2-^2^ ^'-^ - 1' 

Iq, we can obtain the following equations by taking 



(4) 
(5) 



XiX[ — X2X2, 
^iP'dn^l = ^2 PL 2 ^'2^ 



j = 2,...,c. 



Since F has no fully aliased columns, from equation (4), we know that 
and must be the same column of F, thus Xi = X2, and XiPl^ _^X[ = 
XiP^^X'i for J = 2, . . . ,c. Also, since Xi is a column full rank matrix, we 
have P!] ^ — P'd > and thus dji = dj2, for j = 1, . . . , c, that is, d^ = d^. So d^ 



and must be the same row of D since D has no identical rows. Therefore, 
®A and d^ (Ba P are the same column of D' ®a which contradicts 
the fact that they are two different columns of D' (Ba Hence, D' (Ba F as 
well as F ®a have no fully aliased columns. □ 



Proof of Theorem 3. We only prove that there are no fully aliased 
columns between Fi (Bai D' and © F2. (The others can be proved easily.) 
Suppose (Bai and 0^ © are columns of Fi (Bai D' and 0„j © F2, 
respectively, where and d^ = {0,d2i, ■ ■ ■ jdn^i)' are columns of Fi, 

F2 and D' , respectively. Let X and Y be the induced matrices of and 
respectively. Then the induced matrices of d^ (Bai and © are 
[X' , (XPdjj)', . . . , {XPd„ j)'] and Y ® 1^, respectively. From the definition 
of an induced matrix, it is easy to see that XX' ^ yolnilni; where yo is the 
(l,l)th entry of YY' and denotes the ni x 1 vector with all elements 
unity. Thus, from Theorem 1, ®A\ and © 0^^ are not fully aliased. 
Therefore, ©^^^ d^ and 0^ © are not fully aliased. □ 



The following lemma will be used in the proof of Theorem 6. 
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Lemma 3 [Fang, Lin and Liu (2003)]. Suppose is the jth column of 
an F{n, qi - ■ ■ Qm) with induced matrix Xj, j = 1, . . . ,m. Then 

hoB{f\n=triXiX,X'^X,) . 

QiQj 

Proof of Theorem 6. (a) From Theorem 1 and Lemma 3, the induced 
matrices of hi and /12 are Xi 1^X2 and X3 (8^4, respectively. Then we have 



/nod(/ii, /12) = tr[{X[X3X^Xi) ® {X!,Xa'4X2)] 



2^2 



l"'2 



2„2 



l"-2 



triX[XsX'^Xi) tr{X!,XiX'^X2) 

ni=i Qi 

/nod(/^,/^)/nod(/^, /^) 



H —f'NODif^jf^) H ^/nod(/^,/'*)- 

Q2Q4 Q1Q3 

(b) The induced matrices of ®Ai P and are {Xi®X2)P and 

(X3®X4)g, respectively, where P = (P^, P^^_i)', Q = (QJ,, Q^a-i)', 
Pi and Qj are permutation matrices defined by (0, . . . ,gi — 1)P/ = i +a-^ 
(0, . . . , gi - 1) and (0, . . . , 93 - l)Q'j = j +^3 (0, . . . , - 1), respectively, i = 
0,..., gi-l,j = 0,...,<73-l. Let r = P'(X(X30X^X4)Q-ffV4. Then 
from Lemma 3, /nod(/^ ©^3 /^) equals the sum of squares of the 

elements of T. Let W = {wij) = X^Xs, B = {bij) = X'^X^ - n2/(gi(?3)lqi I'q, , 
and note that E'Li Eti = ^i, Ef=i Eti = 0- Then 

93 91 

j=i i=i 

and the (s,t)th entry of T can be expressed as Ylf=iYli=i'^ij^Sitj-, where 
(si, . . . , SgJ and (ti, . . . , ^^3) are some permutations of (1, ... , qi) and (1, . . . , 
(73), respectively. Then 

(93 91 \ 2 / 93 91 \ / 93 91 \ 

i=i i=i / \j=i i=i / \i=i i=i / 

and thus 

/ 93 91 \ / 93 91 \ 
/nod(/' ©^1 f\f ©A3 /') < QiQz E E 4 E E 



,2 

j=l j=l / \j=l i=l 



^193 



Tl?' 

/nod(/^,/^) H ^ 

9193 



/nod(/^,/^'' 
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where the equahty holds if and only if there exist ci and C2 with |ci | + |c2| > 
such that ciWij = 0265,*^ for i = 1, . . . , gi and j = 1, . . . , 53. This means that 
ci Ef=i Eti ^ii = C2 Ef=i Eti bsa^ = 0, and so ci = 0, C2 ^ and 6,^ = 
for ? = 1, . . . , qi and j = 1, . . . , 53. Thus, is orthogonal to 

On the other hand, if is orthogonal to /nod(/^,/^) =0 and 

o</NOD(/'e^i 0^3/') 



/nod(/"'^,/^) + 



n 



2 n 



9193 



/nod(/',/") = 0. 



then the equality holds. 

Similarly, we can obtain that 

/nOd(/' /^ f 0^3 /') = /nOd(/' 0^1 f\ f ®As f) 



< 92Q4 



/nod(/^,/^) + 



9294 



/nod(/^,/^ 



and the equality holds if and only if is orthogonal to f^. Hence, we have 
the assertion. 

(c) The induced matrices of ©^^ and /12 are {Xi ®X2)P and X3 0X4, 
respectively. Let = P' {X[X^® X'^X^) 



712 



^^g^-lg^lg^ and = (wjj) = X(X3, and note that YllLi^t] 
1, . . . , ^3. Then K 



, ^93 ) 5 



where = YlLi'^ijPj-iG- Note that 
/nod(/^ 0.4i /i2) is equal to the sum of squares of the elements of K, the 
(s,i)th entry of Aj is YT=i'^i39sit and {YT=i'^i39sitf < E?ii ^ij EsLi ffst, 
where (si, . . . , SgJ is a permutation of (1, . . . , gi). Then similar to the proof 
in (b), we get 



93 94 91 



91 



91 



/NOD(/^0^./^/.2)<EEE E-^E 

j=l t=l s=l \j=l k=l 
93 91 94 91 

=9iEE^iiEE5fc« 



9kt 



91 



/nod(/^, /^) + 



9193 



/nod(/^,/'' 



where the equality holds if and only if is orthogonal to and 



/nod(/^ 0A /^/i2) <9i 



/nod(/^,/^) + 



n 



2 n 



9194 



/NOD(/^ /^), 



where the equality holds if and only if is orthogonal to Thus, we 
complete the proof of (c) . □ 
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APPENDIX B: SOME SELECTED OPTIMAL SUPERSATURATED 

DESIGNS 

Table 5 

Equidistant designs used in Tables 6-9 
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Georgiou and Koukouvinos (2006) 
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Fang, Ge and Liu (2004) 


16 


15fc (fc = 2,...,6) 


2 


Liu and Zliang (2000) 


16 


15fc (fc = 7,8,9) 


2 


Eskridge et al. (2004) 


16 


5A; (fc = 2,...,6) 


4 
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Georgiou and Koukouvinos (2006) 


20 


19fc (fc = 2,3) 


2 


Liu and Zhang (2000) 


20 


19 


4 


Lu et al. (2002) 


20 


19 


5 


Lu et al. (2002) 


22 


42 


2 


Liu and Zhang (2000) 


24 


46 


2 


Liu and Zhang (2000) 


24 


23 


4 


Lu et al. (2002) 


24 


23 


6 


Lu, Hu and Zheng (2003) 


25 


6fc (/s = 2,...,25) 


5 


Georgiou, Koukouvinos and Mantas (2006) 
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Table 6 

Some selected E{fj^oB)-optimal SSDs constructed by Theorem 3 





mi 




n2 


1X12 


12 


r 


Final resulting SSD^ 


A 


4 


3 


2 


6 


5k 


3 


4A: 


F(24, 2 3 ) 


13fc, A: = 1,2,3 


4 


3 


2 


8 


7k 


4 


6A: 


F(32,2^^ A' ) 


19fc, fc = 1, . . . ,6 


4 


3 


2 


9 


4:k 


3 


3k 


F{36,2^^*'3'^^) 


lOfc, k = 2t,t = l...,6 


6 


5 


3 


6 


10 


2 


3 


F(36,2"^3 ) 


19 


4 


3 


2 


10 


9k 


5 


8k 


_F(40, 2 5 ) 


25fc, fc = 1, . . . ,4 


4 


3 


2 


12 


Ilk 


3 


8k 


F(48,2 " 3 ) 


27fc, fc = 1, . . . ,5 


6 


5 


3 


8 


7k 


2 


2k 


F(48,2^''3'^" ) 


IS/bj k — 2j . . . ^ 5 


6 


10 


2 


8 


7k 


4 


3k 


^(48,2^^" 4' ) 


31fc, fc = 2,4,6 


4 


3 


2 


12 


Ilk 


6 


10k 


F(48,2'^°*6^^'°) 


31fc, fc = 1,2,3 


6 


5 


3 


8 


7k 


4 


3k 


F(48,3^ M^*") 


16fc, fc = 1, . . . , 6 


6 


10 


3 


8 


7k 


2 


k 


F(48,2'^''3^°'=) 


13fc, fc = 3,4,5 


6 


10 


3 


8 


14fc 


4 


3k 


F(48,3^" 4^*") 


32fc, fc = 1,2,3 


6 


10 


2 


9 


16 


3 


6 


F(54,2^^°3^'*) 


64 


4 


3 


2 


14 


13fc 


7 


12fc 


F(56,2^^''7"'°) 


37fc, fc= 1,2 


4 


3 


2 


15 


28 


3 


20 


F(60,2^''"3 **) 


68 


6 


5 


3 


10 


18fc 


2 


5k 


F(60,2^**''3^^ ) 


33/i(5 — 2^ 3 


4 


3 


2 


15 


7k 


5 


6k 


F(60,2^'^''5^'') 


19fc, fc = 2,...,13 


6 


10 


2 


10 


9k 


5 


Ak 


F(60,2**°'-5^'-) 


41fc, A: = 2,3,4 


6 


5 


3 


10 


9k 


5 


Ak 


F(60,3'^°'=5^*) 


21fc, fc = 1,...,4 


10 


9 


5 


6 


5k 


3 


k 


F(60,3^'=5'"'*) 


lOfc, fc = 2,3 


6 


10 


3 


10 


18 


5 


4 


F(60,3^^°5^®) 


42 


4 


3 


2 


16 


5k 


4 


4fc 


F(64,22**=4^*=) 


13fc, fc = 2,...,16 


4 


3 


2 


18 


17k 


3 


12fc 


. . , »70 I-. 1 *7i[-. y 

F(72,2'^2*3"*) 


41fc, fc = l,2 


6 


10 


2 


12 


Ilk 


3 


Ak 


F(72,2*" 3 ) 


43fc, fc = 2, ...,5 


4 


3 


2 


18 


34 


6 


30 


F(72,2^*°6^*) 


94 


6 


10 


2 


12 


22 


6 


10 


F(72,2"'""6 ) 


102 


9 


4 


3 


8 


7k 


4 


6k 


T— • / -^70i* .7i*\ 

F(72,3^^''4^'') 


25fc, fc= 1,...,6 


U 


c; 


Q 
O 


1 9 




« 

u 




^^79 Q75fecllfe\ 


ZUA,, A, O 


10 


18 


2 


6 


5k 


3 


2k 


F(80,2^2'=3^'=) 


37fc, fc = 2,3 


10 


18 


2 


8 


7k 


4 


3k 


F(80,2i°*'=4^'=) 


55fc, A: = 2, 4, 6 


4 


3 


2 


20 


19 


5 


16 


F(80,2^^5'^) 


51 


10 


9 


5 


8 


Uk 


4 


3k 


F(80,4^*'=5"^'=) 


29fc, fc = l,2,3 


6 


10 


2 


14 


26 


7 


12 


F(84,22«72S) 


122 


6 


5 


3 


14 


26 


7 


12 


F(84,3i«°72'^) 


62 


6 


10 


2 


15 


7k 


5 


3k 


F(90,2'^°'=5'''') 


31fc, k = 2t,t = 2,...,6 


6 


5 


3 


15 


7k 


5 


3k 


F(90,3'"^'=5'^'') 


16k, k = 2...,13 


9 


4 


3 


10 


9k 


5 


8k 


F(90,3^®*5^*) 


33k, fc = 1,...,4 


6 


5 


3 


16 


15k 


2 


Ak 


F(96,2i'5fe360fc) 


27fc, fc = 2,...,9 


4 


3 


2 


24 


23 


4 


18 


F(96,2"«423) 


59 


6 


10 


2 


16 


5k 


4 


2k 


F(96,2''°'=4'^'=) 


21fc, fc = 4,...,16 


4 


3 


2 


24 


23 


6 


20 


F(96,2i20g23) 


63 


6 


5 


3 


16 


5k 


4 


2k 


F(96,3^°'=4'5'=) 


llfc, fc = 3,...,16 


4 


3 


2 


25 


6k 


5 


5k 


F(100,2^'"=5'"=) 


16A;, A; = 2t,t = 2, ...,12 



A is the constant coincidence number of the final resulting SSD. 
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Table 7 

Some selected -optimal SSDs constructed by Theorem 3 





mi 




7l2 


771-2 


Q2 




Final rpsiiltiriD- SSD^ 




4 


3 


2 


6 


5k 


3 


6fc 


^(24,2^^*3'^*) 


39k, k^ 1,2,3 


4 


3 


2 


8 


7k 


4 


12fc 


F(32, 2^2*4^*) 


76k, k=l,...,6 


4 


3 


2 


9 


16k 


3 


18fc 


F(36,2io»*3i'^'=) 


120fc, fc=l,2,3 


6 


10 


2 


6 


10 


3 


6 


^^(36,2120310^ 


126 


6 


5 


3 


6 


10 


2 


2 


F(36,2i°3^°) 


38 


4 


3 


2 


10 


9k 


5 


20fc 


F(40,2i2°'=5^'=) 


125fc, fc = l,...,4 


4 


3 


2 


12 


Ilk 


3 


12fc 


F(48, 2^2^311;=-) 


81fc, fc=l,...,5 


8 


14 


2 


6 


10 


3 


6 


F(48,2i«*3'") 


174 


4 


3 


2 


12 


life 


6 


30fc 


F(48,2i»°*6"'=) 


186fc, fc=l,2,3 


6 


10 


2 


8 


7k 


4 


6k 


F(48,2i2°'=4^'=) 


124fc, fc = l,...,6 


6 


5 


3 


8 


7k 


4 


4k 


f (48,3'*'"= 4^*=) 


64fc, fc= 1,...,6 


6 


10 


3 


8 


7k 


4 


2k 


^(4g^ 360*47* j 


64fc, fc = 2, ...,6 


6 


15 


3 


8 


42 


4 


8 


^(433360442-) 


384 


8 


7 


2 


6 


5k 


3 


6k 


F(48,2'**'=3'5'=) 


87fc, fc = l,2,3 


8 


7 


4 


6 


5k 


3 


k 


F(48,42**'=3'"=) 


31fc, fe = 2,3 


4 


3 


2 


14 


13 


7 


42 


F(56,22'^27i^) 


259 


4 


3 


2 


15 


28 


3 


30 


F(6Q,2i«''32») 


204 


4 


3 


2 


15 


14 


5 


30 


F{60,2^''°5'^^) 


190 


6 


10 


2 


10 


9k 


5 


lOfc 


F(60,22'''"'5»'=) 


205fc, fc = l,...,4 


10 


18 


2 


6 


10 


3 


6 


F(6Q,22i63i'') 


222 


4 


3 


2 


16 


5k 


4 


8fc 


(64, 2 4 ) 


52k, k = 2, ... ,1b 


8 


14 


2 


8 


7k 


4 


6k 


^(64,2^'***'= 4^'') 


172k, fc = l,...,6 


4 


3 


2 


18 


17 k 


3 


18k 


F(72,2i''**3i'^*) 


123fc, fc = l,2 


8 


7 


2 


9 


8k 


3 


9k 


F(72,2i2'^'"3*'=) 


132fc, fc = 2,4,6 


9 


8 


3 


8 


7k 


4 


4k 


F(72,3'""=4'"=) 


lOOfc, fc = l,...,6 


9 


16 


3 


8 


7k 


4 


2k 


F(72,3'"*'=4'''=) 


lOOfc, fc = 2,...,6 


6 


5 


3 


12 


Ilk 


6 


10k 


F(72,3i'50fcgiifc) 


156fc, A; =1,2, 3 


8 


7 


2 


10 


9k 


5 


20k 


F(80,228'"=5^'=) 


285fc, fc = l,...,4 


10 


18 


2 


8 


7k 


4 


6k 


F(80,22i6'=4^'=) 


220fc, fc = l,...,6 


6 


5 


3 


14 


13fc 


7 


14k 


F(84, saio^risfc) 


217fc, fc = l,2 


6 


5 


3 


15 


7k 


5 


5k 


F(90,3''5'=2^*) 


80fc, fc = 2,...,ll 


6 


10 


3 


16 


30 


4 


8 


F(96,32'"^43'') 


264 


10 


18 


2 


10 


9k 


5 


10k 


F(100,23'""=2^'=) 


365fc, fc = l,...,4 


9 


4 


3 


12 


Ilk 


2 


4k 


F(108,3^**'=2"'=) 


58fc, fc = l,...,12 


10 


18 


2 


12 


Ilk 


3 


6k 


f (120,22i'*'=3"'=) 


225fc, fc = l,...,5 


8 


7 


4 


14 


26 


7 


14 


F(112, 43^242'') 


406 


8 


14 


2 


16 


5k 


4 


4k 


i^(128,2"2'=4*'') 


116A:, fc = 2, ...,16 


9 


4 


3 


15 


7k 


5 


10k 


F(135,3'2'"=5^'=) 


125fc, fc = l,...,5 



^F{n,n2,qr"''q"^'). 

uj is the constant natural weighted coincidence number of the final resulting SSD. 
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Table 8 

Some selected E^fi^ou) -optimal SSDs constructed by Theorem 5 



m 


TTli 


Qi 


n2 


m2 


Q2 




Q3 


m4 


<74 




r4 


Final resulting SSD^ 


A 




4 


3 


2 


6 


10 


2 


3 


2 


5 


3 


6 


4 


F(24,2^'^3^''4^°) 


42 




4 


3 


2 


6 


10 


3 


3 


2 


5 


3 


8 


2 


F(24,2***3^''6^'') 


36 




4 


3 


2 


6 


15 


3 


3 


2 


5 


3 


12 


3 


f (24, 2^2345545-) 


54 




4 


3 


2 


8 


7k 


2 


3 


2 


7 


4 


4fc 


2k 


f (32, 224*4^^'=) 


29fc, k = l,.. 


..,5 


4 


3 


2 


8 


21 


2 


3 


2 


14 


4 


12 


3 


F(32,2''^4^^^) 


87 




4 


3 


2 


8 


28 


2 


3 


2 


7 


4 


16 


8 


F(32,29«4308) 


116 




4 


3 


2 


8 


21 


4 


3 


2 


7 


4 


18 


2 


F(32,2i°*^45«8'^^) 


71 




6 


10 


2 


6 


10 


2 


10 


2 


5 


3 


12 


12 


J?(36,22«3i«°4i00) 


196 




4 


3 


2 


9 


8k 


3 


3 


2 


4 


3 


6k 


4fc 


F(36,236'=3'*8'=62'*'=) 


36fc, k — 1, . 


..,6 


6 


10 


2 


6 


5k 


3 


10 


2 


5 


3 


8k 


3fc 


f(36,2i''°'=3^5*g50fe-) 


99k, fc= 1,: 


2,3 


6 


10 


2 


6 


10 


3 


10 


3 


10 


2 


8 


6 


F(36,2i203240gi00) 


148 




6 


5 


3 


6 


10 


2 


5 


3 


10 


2 


3 


8 


F(36,2i'5°3*^6'^°) 


99 




6 


5 


3 


6 


5k 


3 


10 


2 


5 


3 


2fc 


2fc 


F(36,2''°'=33°*=925'') 


31fc, fc = 2,3 


6 


5 


3 


6 


15 


3 


10 


3 


10 


2 


3 


6 


F(36,2i203^°9^5) 


93 




6 


10 


2 


8 


7 


2 


10 


3 


7 


2 


4 


18 


^(48,225231204™) 


178 




6 


5 


3 


8 


14 


2 


5 


3 


14 


2 


4 


12 


^(4g 2336360g70) 


194 




6 


5 


3 


8 


21 


2 


10 


3 


14 


2 


3 


18 


j,(4g 2504390gl05) 


291 




6 


5 


3 


8 


28 


2 


5 


3 


21 


2 


8 


16 


F(48,28^23i20gi40) 


388 




6 


5 


3 


8 


7k 


4 


5 


3 


7 


2 


3fc 


4fc 




44fc, k = l,.. 


..,6 


6 


10 


2 


8 


7k 


2 


10 


2 


14 


4 


8fc 


3fc 


f (48,2^^°''423**) 


134fc, fc=l,. 


..,5 


6 


5 


3 


8 


7k 


2 


10 


2 


7 


4 


2k 


4fc 


f (48,2'""=4"2'=6^'"=) 


51fc, fc = 2,. 


..,5 


6 


5k 


3 


8 


21 


2 


10 


2 


7k 


4 


6k 


12 


F(48 2-'-2("'4336*gl05*'j 


153fc, k = l, 


2,3 


6 


10 


2 


8 


7k 


2 


5 


3 


7 


4 


8k 


6fc 


f (48,3^20*4238*-) 


94k, k^l,. 


..,5 


6 


5 


3 


8 


14 


2 


5 


3 


7 


4 


4 


8 


^(43^360 4224g70) 


82 




6 


5k 


3 


8 


21 


2 


5k 


3 


7k 


4 


6 


12 


F(48 390'=4336fegl05fe-j 


123fc, k = l. 


2,3 


6 


10 


3 


9 


8 


3 


10 


2 


8 


3 


6 


8 


^(54,21203192980) 


128 




10 


18 


5 


6 


10 


2 


18 


2 


5 


3 


6 


32 


F(60, 2216348010180) 


276 




6 


10 


3 


10 


9 


5 


5 


3 


18 


2 


8 


4 


F(60,2i**3i20l5''0) 


114 




6 


10 


3 


10 


9 


5 


5 


3 


9 


5 


8 


2 


^(60,31205901590) 


60 




8 


7 


2 


8 


7k 


2 


7 


2 


7 


4 


12fc 


4fc 


F(64, 2i'^*''4i''i'') 


121fc, k=l,. 


..,5 


8 


7k 


2 


8 


14 


2 


14 


2 


7 


4 


12fc 


8fc 


7?(64|2^^'"=4^22fc) 


242k, k=l,. 


..,5 


8 


7 


2 


8 


7k 


4 


21 


4 


7 


2 


2k 


4fc 


_F(64,256'°4i^*''8'i'"') 


73k, k^l,. 


..,6 


8 


7k 


2 


8 


14 


4 


21 


4 


14 


2 


4fc 


4fc 


F(64 2II2'' 4336*898*:) 


U6k, fc=l,. 


..,5 


8 


7k 


2 


8 


21 


4 


21 


2 


7 


4 


18fc 4fc 


f (64,2''''6'=4ii2*8i'''"') 


415fc, k=l,. 


..,5 


8 


7k 


2 


9 


16 


3 


7k 


4 


16 


3 


12 


4fc 


-F(72 3l92*4336*gll2*) 


160fc, k=l,. 


..,5 


10 


9k 


5 


8 


14 


2 


9k 


5 


7fc 


4 


2 


16 


f(80,'4*4»'=5^»'=10i2«'=) 


136k, k=l,. 


..,4 


10 


9k 


5 


8 


21 


2 


9k 


5 


14fc 


4 


3 


12 F(80,4672fc5i35fe^oi««'=) 


204fc, k = l. 


2,3 


10 


9k 


5 


8 


14 


4 


9k 


5 


7k 


2 


3 


16 F(80,2224'=5i35fc20i26fe) 


Ulk, k=l,. 


..,4 



t F (nma , (gi 92 ) 1 '"^ Qa""'' ) ; "1 = ^3 , na = n4 . 

A is the constant coincidence number of the final resulting SSD. 
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Table 9 

Some selected -optimal SSDs constructed by Theorem 5 



m 


mi 


<7i 


^2 


m-2 


92 




Qs 


m4 


Q4 




r4 


Final resulting SSD^ 






4 


3 


2 


5 


5 


3 


3 


2 


5 


3 


12 


2 


pz-r,^ r,72o30(:.15\ 

-1. I Zi^:, Zi O U f 


108 




4 


Q 
O 


2 


e 

o 


I n. 


2 


Q 

o 


2 


7 


4 


8A; 


2k 


p/or, y48k,77k\ 


1 1 (^i, Z:- — 1 


. . , 


Q 


8fc 


Q 
O 


4 


Q 
O 


2 


8fc 


Q 
O 


Q 
O 


2 


4 


18k 


p/of? 9l08fcq96fc/?24fc\ 


Z±IJr»(, ft. — _L , . 


. . , D 


u 




Q 
O 


R 
U 


1 n 


Q 
O 


1 n 


9 
z 


1 n 


Q 

o 


1 Si- 


Un- 


j^da o360fcQl80fcn50fc\ 
-f i^ou, z o y J 


— 1 

OOOft. , n- — -L 5 


9 

Z, ,3 







9 


g 
o 


7 
i 


9 


1 n 

lU 


9 
z 


1 /I 


/I 


1 f\ 
10 


Q 
O 




OoO 










Q 
O 


g 
o 


7 
t 


9 


1 n 


9 
z 


1 /I 


/I 


D 


Q 
O 


p/.o 91204168^^35^, 

-T i^'^o, z D J 


oUO 




« 
u 


o 


Q 
O 


g 

o 


7 


9 




Q 
O 


1 4 


9 
z 


A 


1 8 


F(AR 9504o60^35n 
r ^^^o, z \j ) 


Ooz 




« 
u 




Q 
O 


g 

o 


7Jc 

t 


9 
z 


c; 
o 


Q 
O 






Ah 


Q 
O 


r?{ AO Q60fc/il68fcf>35fc\ 




9 ^ 
z, a 


« 
u 




Q 
O 


g 

o 


1 A 


9 
z 


1 n 


9 
z 


7 




1 9 
±z 


1 9 
±z 


PC48 92404336f.70N 


612 




R 

u 


K 
tj 


Q 
O 


g 

o 


14 


2 


o 


Q 
O 


7 


4 


g 

o 


12 


PC48 31204336570-, 


492 




« 
u 


K 


Q 
O 


g 

o 


21 


2 


1 n 


2 


42 


4 


1 S 


Q 
O 


F(A8 9360 4504 105 ^ 


918 




g 


5 


3 


8 


28 


2 


10 


3 


42 


4 


8 


4 


Fi48 32«4672gl40-, 


984 




6 


5 


3 


8 


7 


4 


10 


2 


7 


4 


18 


4 


^(4g 23604112^235) 


484 




6 


5 


3 


8 


7 


4 


5 


3 


14 


2 


12 


12 


^(4g|23363l80;^235) 


528 




6 


5 


3 


8 


7 


4 


5 


3 


7 


4 


12 


4 


^(4g^3i804ii2^235) 


304 




6 


5 


3 


8 


14 


4 


10 


3 


21 


2 


12 


16 


^(4g 26723860^2™) 


1056 




6 


5 


3 


8 


7fc 


4 


10 


3 


7 


4 


6k 


4k 


-F(48 3180^4112*: -j^235'=) 


304fc, k = l,. 


..,6 


6 


5 


3 


8 


21 


4 


10 


3 


7 


4 


18 


12 


^(4g^35404336^2l05) 


912 




6 


10 


2 


9 


4 


3 


10 


2 


4 


3 


18 


12 


^(54 23603l44g40) 


528 




10 


9fc 


5 


6 


5 


3 


9fc 


5 


5k 


3 


3 


20 


F(60, 3300*5135* IS''^*) 


450fc, k=l, 


2,3 


8 


7 


2 


8 


7 


2 


14 


2 


7 


4 


12 


4 


F(64,23364i6i) 


484 




8 


7 


2 


8 


14 


2 


21 


2 


7 


4 


16 


8 


F(64, 2^^24322-, 


968 




8 


7 


2 


8 


21 


2 


28 


2 


7 


4 


18 


12 


F(64, 210084483) 


1452 




8 


7 


2 


8 


7k 


4 


7k 


4 


7k 


2 


12 


16 


_F'(64 222'*'=4336feg49fe-j 


584fc, fc = l,. 


..,5 


8 


7 


2 


8 


21 


4 


21 


4 


28 


2 


12 


12 


F(64,2'5^2 4ioos8i47) 


1752 




8 


7 


2 


9 


8 


3 


21 


2 


8 


3 


18 


8 


f(72,2^563l92g56) 


984 




8 


7 


2 


9 


8fc 


3 


21fc 


4 


16 


3 


3 


4k 




480fc, fc = l,2 


8 


7 


4 


9 


4 


3 


7 


2 


8 


3 


18 


12 


F(72, 2^523288^228) 


552 




8 


7 


4 


9 


4fc 


3 


7k 


4 


8k 


3 


3 


12 


F{72 3288*1484*: ■j^228fc'j 


384fc, k = l,. 


..,6 


9 


8fc 


3 


8 


7 


2 


8k 


3 


21 


2 


8 


18fc 


F{72 2^56fc3l92fcg56fc^ 


984k, k^l,. 


..,6 


9 


8fc 


3 


8 


7 


2 


8k 


3 


21 


4 


8 


3k 


F{72 3l92fe4252feg56fe'j 


480fc, fc = l,. 


..,6 


9 


8fc 


3 


8 


7 


4 


16 


3 


21 


2 


12fc 12fc f^(72, 2^^**3'''^^* 12^5^*) 


1104fc, fc=l, 


...,4 


9 


8fc 


3 


8 


7 


4 


16 


3 


14 


4 


12fc 


3k 


F{72 35'''6*i4l68*! 2^256*:^ 


768fc, fc = l,. 


..,4 


10 


9 


5 


8 


7k 


2 


18 


2 


7k 


4 


lOfc 


20 


-F(80 23^0*^4560/2 "1^063*;^ 


950fc, fc = l,. 


..,5 


8 


14 


2 


10 


9k 


5 


28 


2 


9k 


5 


12fc 


2 


F(8o!2«^2;=59o;=;LOi26fe-) 


360fc, fc = 1,. 


..,4 


10 


9fc 


5 


8 


7 


4 


9k 


5 


28 


2 


6 


20fc F(80,2"2»'=52^°'=20«3'=) 


1410fc, k=l, 


...,5 


8 


7fc 


2 


10 


9 


5 


7k 


4 


9 


5 


20 


2k 


^(gQ^4560fc590fc^063fe) 


680fc, fc = 1,. 


..,5 



^ F(ni na , (gi ga ) "1 '"^ Qa"'''' ) ; '^i = ^^3 , na = n4 . 

a; is the constant natural weighted coincidence number of the final resulting SSD. 
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